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1. Introduction 

In this paper we investigate liow symmetry can be broken in some variational problems 
when an external potential or a weight enters in competition with a nonlinearity. An 
archetypal example for such issues is the question of symmetry of optimal functions in 
Caffarelli-Kohn-Nirenberg inequalities. While all terms are invariant under rotation 
around the origin, it is known that optimizers are not always radially symmetric. 
Caffarelli-Kohn-Nirenberg inequalities, also known as Hardy-Sobolev inequalities, is a 
particularly simple setting for the study of symmetry breaking because weights and 
nonlinear terms have simple homogeneity properties, so that Euler-Lagrange equations 
inherit scaling properties that allow to further simplify the symmetry issues. Still, ranges 
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of parameters for which optimizers are radially symmetric have not been completely 
determined yet. 

Symmetry breaking issues are present in many areas of physics involving partial 
differential equations: quantum mechanics, mean field models, equations for phase 
transition, ferromagnetism, mechanics, etc. Symmetry breaking occurs when antagonist 
effects are competing, like weights or potentials (or coupling with other fields) on the 
one hand and nonlinearites on the other hand. Various mathematical methods are 
available either for proving symmetry (uniqueness, comparison techniques based for 
instance on moving planes methods, symmetrization: see for instance [H [21 E]) or for 
proving symmetry breaking (multiplicity and bifurcation, energy, spectral methods). 
However, threshold cases are not characterized even in the simplest cases. 

A simple mechanism which can break symmetry is the instability of the symmetric 
extremals, that is, the case where the extremals among radially symmetric functions are 
not local minima in the larger space of functions with no symmetry assumption. In 
the case of the Caffarelli-Kohn-Nirenberg inequalities, this instability has been studied 
in several papers (see [H El E]) and the corresponding region of symmetry breaking is 
delimited by a curve that we will call for brevity the curve of Felli- Schneider. However, 
it has been proved in [3] that symmetry breaking can occur even in a range of parameters 
for which the symmetric extremals are stable, that is, in cases where they are strict local 
minima. In order to understand this phenomenon, and symmetry breaking in general, 
we study the solution set associated to the Euler-Lagrange equations corresponding to 
a minimization problem associated to the Caffarelli-Kohn-Nirenberg inequalities. For 
those equations, we investigate the bifurcation of non-radially symmetric solutions from 
radially symmetric ones. The two theoretical contributions of the present paper are 
an asymptotic analysis of the branches for large values of the bifurcation parameter, in 
Section |2| and a detailed expansion of the non-radial solutions in a neighborhood of 
the bifurcation point on the branch of radial extremals, in Section |3| Both results are 
consistent with known numerical results and give a significant insight about the local 
behaviour of the solutions, either around the bifurcation point or asymptotically, that 
is numerically difficult to compute by direct methods. Let us give some further details. 

We shall consider a subfamily of the Caffarelli-Kohn-Nirenberg inequalities which, 
for a given dimension d > 3, depend on two exponents, p G (2, 2*] with 2* := 2 d/{d — 2) 
and 9 G ['(9(p, d), 1] with 'd{j),d) := d {p — 2)/{2p), and on a parameter A > (see Section 



2.1). We shall also consider the two-dimensional case: d = 2, p G (2,2*) with 2* = oo. 
Partial symmetry results are known from [71 [HI E] and it has also been shown in fiUl [3] 
that, in the set of parameters (p. A), a curve separates the region of symmetry from 
the region of symmetry breaking. A priori estimates can be used to localize this curve 
(primarily see [11], and references therein), but the curve itself has not been properly 
characterized yet. Our main purpose is to study branches of solutions to the Euler- 



Lagrange equations when 6 = 1 (see Section 2. 1| below) and first check the conjecture 
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that the threshold is exactly given by the curve of Felli-Schneider, as was strongly 
suggested by [11] . We also aim at understanding the symmetry breaking observed in |3l 
Theorem 5] when 9 < 1. This can be done by reparametrizing the solutions set found 
in the case 9 = 1. 

Numerical branches of solutions have been computed in great detail in [12]. The 
scheme allows us to construct branches of non-symmetric solutions which are candidates 
to be minimizers. Even if it cannot be granted that the computed solutions are 
extremals, their qualitative properties fully account for all observations, including the 
instability results in the range of parameters for which the symmetric extremals are 
stable. Based on these computations, two global scenarii can be distinguished at 
a qualitative level: for a given dimension d > 2 and a fixed exponent p, either the 
branch of non-symmetric solutions is monotone increasing for all admissible values of 
the parameter 9, or not. Numerically, we observe that the first scenario always holds 
for 9 close to 1 while the second one is observed if p is not too large and only when 9 is 
small enough. This picture is reinforced by two observations: 

(1) What happens at the bifurcation point is difficult to characterize numerically using 
the scheme of [12] • However, the expansion of Section [s] shows that the branch is, 
locally monotone increasing if and only if 9 is bigger than some threshold value. 

(2) In the limit case 9 = ^{p^d), the branch of non-symmetric solutions has an 
asymptote. As in [13], the limiting energy can be compared with the energy of 
symmetric solutions at the bifurcation point. In case it is higher, we numerically 
observe that the branch is monotone increasing for all values of 6* G ['(9(p, c?), 1], but 
not in the opposite case. 

These two observations are consistent with the global qualitative scenarii. The numerical 
agreement of the threshold value for (1), which is local and restricted to a neighborhood 
of the bifurcation point, and the threshold value for (2), which is purely asymptotic, is 
striking. From a numerical point of view, this is our main result. This suggests that 
there is no other scenario than the two described above. 

Our paper is organized as follows. Some known facts and preliminary estimates 
are collected in Section [2| including the method of the reparametrization of the branch 
of solutions for 6' < 1 by the solution set corresponding to ^ = 1. The analysis of the 



asymptotic regime as A — )■ oo is performed in Section 2.4 when 9 G ['i9(p, d), 1] and 
the corresponding result. Proposition [11 is new when 6* 7^ ■i9(p, d) or 1. The expansion 
at the bifurcation point is done in Section [3} at a formal level. Indeed we make an 
ansatz for the non-symmetric solutions, and study the behavior of the branches near 
the bifurcation point within the corresponding class of functions. Most likely the ansatz 
can be justified to the price of complicated but elementary a priori estimates. This 
is the core of our contribution and it goes well beyond the stability analysis that can 
be found in [U E] E] . Numerical results have been collected in Section 111 and discussed 
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there. Section [5] is devoted to some concluding remarks. 

2. Preliminaries and theoretical results 

2.1. Caffarelli-Kohn-Nirenherg inequalities 



For any dimension c? G N*, let us consider the set V of all smooth functions which are 
compactly supported in W^. Define the numbers 



p — 2 d — 2 

d{p, d) := d — — , Oc := -—— , A(a) := (a - acf , p(a, h) : = 



2d 



d-2 + 2{b-a)' 



We also set 2* := -^^ if rf > 3 and 2* := oo ii d = 1 or 2. For any a < Oc, we consider 
the following Caffarelli-Kohn-Nirenherg inequalities, which have been introduced in [H] 
(also see 0): 

Let b e {a + 1/2, a + 1] and 9 E (1/2, 1] if d = 1, b E (a, a + 1] if d = 2 and 
b E [a, a + 1] if d > 3. Assume that p = p{a, b), and 9 E [■(9(p, rf), 1] if d > 2. Then, there 
exists a finite positive constant Kckn(^, A,p) such that, for any w eT), 



I i-fe ||2 ^KcKN(^,A(a),p) 20 



I |-(a+l) ||2(l-6») 

l-^l '^|lL2(IRd) 



(1) 



We denote by Kcj^j^(6', A(a),p) the best constant among all radial functions. We 
recall that this constant is explicit (see [6l Lemma 3]) and equal to 



with K 



CKNl'^5 



l.P)- 



(p-2)^ 



(2e-i)p+2 



P-2 
2p 



{2e-i)p+2 


e 


4 


6-p 
2p 


2pe 




p+2 













v^r(^; 



£-2 



According to |5], the Caffarelli-Kohn-Nirenberg inequalities on W^ can be rewritten 
in cylindrical variables using the Emden-Fowler transformation 



s = log Ixl , u = —- E S'^~^ , u(s, u) 
\x\ 



Ixl"" "wix) 



and are then equivalent to Gagliardo-Nirenberg-Sobolev inequalities on the cylinder 

C := M X S"^"-^ that can be written as 



\u\ 



LP(C) 



< Kckn(^,A,p) (\\Vu\\l2^c)+^Mh 



(c) 



\u 



|2(i-e) 
Il2(c) 



(2) 



for any u E H^{C). Here we adopt the convention that the measure on S'^ ^ is the 
uniform probability measure. 
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Radial symmetry of some optimal functions in ([I]), or symmetry to make it 
short, means that there is an optimal function in ^ which only depends on s and 
Kckn(6',A,p) = Kq-^^{6,K,p). We say that there is symmetry breaking if and only 
if Kckn(6',A,p) > K^j^i^(6', A,p). Notice that on the cylinder the symmetric case is 
equivalent to the one-dimensional Gagliardo-Nirenberg-Sobolev inequality 






for any u G iJ^(]R). Symmetry breaking of course makes sense only if rf > 2 and we will 
assume it is the case from now on. Let us summarize known results. Let 

^ , ^, d-l{2e-l)p + 2 ^ ,,. (rf-l)(6-p) 
AfsP,^ :=4^— ^ -^ and A, p := ^— y^^— ^ . 

Symmetry breaking occurs for any A > Afs(p, ^) according to |4, iGj (also see for 
previous results and p^ if rf = 2 and ^ = 1). If ^ = 1, from [11], we know that 
symmetry holds for any A < Aj,(p). Moreover, as shown in [T^ (TUl E], there is a 
continuous curve p H- As(]9) with limp^2+ As(p) = oo and As(p) > a1 for any p G (2,2*) 
such that symmetry holds for any A < Ag and there is symmetry breaking if A > Ag. 
Despite the fact that the exact shape of A^ is not known, it can be proved that we 
have limp^2* As(p) = a^ if d > 3 and, if rf = 2, limp_j.oo As(p) = 0, or more precisely, 
limp_^ooP^As(p) = 4. Moreover, we also know from [TJ Theorem 3.1] that symmetry 
holds ifA<rfVp2. 

According to [TB], existence of an optimal function is granted for any 9 G (^(p, rf), 1), 
but only if Kckn(^,A,p) > Kgn(p, c^) when 6 = 'd{p,d), where Kgn(p, c^) is the optimal 
constant in the following Gagliardo-Nirenberg-Sobolev inequality 

|L,||2 ^ Kgn(P,'^) 11^ ||2tf(p,d) II ||2(l-tf(M)) v„ c uim'i^ 

||«|lLP(Rd) ^ |gd-l|(p-2)/p ll^"llL2(Md) ll"llL2(Md) V U t n I^U:^ J . 

A sufficient condition for existence of extremals can be deduced, by comparison with 
symmetric functions, namely K*Qy^-^{'d{j),d),K,p) > Kgn(p, c^), which can be rephrased 
in terms of A as A < Aqp^(p, d) for some non-explicit (but easy to compute numerically) 
function p i— )■ h*Q^{p,d). When 6 = ^{p,d) and A > AQ^{p,d), extremal functions (if 
they exist) cannot be symmetric and in the asymptotic regime p — )■ 2+, this condition is 
weaker than A > Aps{p,'&{p,d)). One can indeed prove that \im.p^2+ -^fs {p,^{p,d)) > 
limp^2+ Aqj,^ (p, d) . Hence, for 9 G {{}{p,d),l], close enough to ^{p,d) and p — 2 > 0, 
small (but numerically not so small, actually, as shown in ^3\ ; also see |3], Section 5] for 



estimates and Section 4.3 for some plots), optimal functions exist and are not symmetric 
if A > Aqj^(p, d), which is again a less restrictive condition than A > Afs(p, 9). See [3] 
for proofs and [TH] for a more detailed overview of known results. 

In this paper we study perturbatively the non-symmetric solutions lying in the 
first branch bifurcating from the branch of symmetric extremals and show that they 
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explain all phenomena of symmetry breaking known or observed so far, including cases 
in which the symmetric extremals are stable. Of course, it is not clear that all extremals 
for Caffarelli-Kohn-Nirenberg inequalities lie in those branches, even if probably that 
is the case. In this paper we provide a complete description of the branch around 
the bifurcation point, based on a formal asymptotic expansion. This clarifies the local 
behavior of the branch and accounts for all phenomena numerically observed in [12] . 



2.2. The case of symmetric extremals 

We start with the symmetric case for 9 = 1 and adapt the computations that can be 
found in [6] (also see [H] and the Appendix). Consider the equation 



-{j>-2fw" + Aw-2p\w\P-'^w = Q in M. (3) 

2 

'unction w{s) := (coshs) p-^ is, up to translations, the unique positive solution 
. As a consequence, the function u{s) = (^p/i) w (^^ ^/^s^ is the unique 



The 
of (|3 

solution of 

- u" + fiu= \u\P-^u in M. (4) 

Now let us turn our attention to the case 6 < 1. Let 
Ic |V«P dy 



t[u] :-- 



Ic w^ dy 



In the cylinder C, solutions of the Euler-Lagrange equation satisfied by the extremals 
of Q satisfy 



- eAu + 



(1 - 6) t[u] +A u = uP-^ (5) 



after a multiplication by an appropriate constant. Such a multiplication has indeed no 
incidence neither on the value of t[u] nor on the value of the quotient 

nei 1 _ Ucl^u^d y + Aj^u^dyY {J^u^ dyY'' 



UcUPdyy 

since u ^-^ t[u] and Q^N] ^^^ both zero-homogeneous in u. Multiplying ([5| by m and 
integrating on C, we find that /^ |VmP dy + A f^u^ dy = f^u^ dy. As a consequence, 
we have 

n^ I 1 _ II ii^p-2 II \\'2{^-d) 
^aFJ — \\'^\\lp{c) II'"IIl2(c) • 

Denote by m* the solution of (|5| depending only on s. The same formulae as for u 
also hold for u^. To relate Kckn(^,A,p) and Kqj^j^{9,A,p) we have to compare Q^N 
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with Q^[m*]. Then, either Kckn(6',A,p) = K*qy^^{6,K,p) and symmetry holds, or the 
inequahty 

is strict and then, there is symmetry breaking. The symmetric optimal function u^ for 
6 < 1 can be explicitly computed. Up to multiplication by a constant, u^ solves 

-eul + riu, = ul-\ 

with Tj = {1 — 6) t[u^] + A. After multiplying the above equation by m*, integrating with 
respect to s e M and dividing by Jj^ u1 ds, we find 

t[u,]+k= f * 



where m*(s) = Aw{Bs), for all s G M, w solves (3), A = {^)''~^ and B 
From this expression, as in pjj, we deduce that 



P-2 /W 

2 V e- 



. J2 p-2r] _^ L^Pds Ip 2pr] 



t\uJ = B'-^ = "^ 4 and '^ : . = AP 



I2 p + 2 9 j^ulds I2 p + 2 



where for all g > 2, 1^ := Jj^ |w(s)|'^cis, and J2 := /jg |ti;'(s)p cis (see Appendix A.l for 
details). This provides the identity 

P-2 t] ^^ ^ 2pr/ 
p + 2 e p + 2 

and uniquely determines i] = 129^1-! ^2 A- As a consequence, we have 

twj = i-- A. 

^ ^ (2^-l)p + 2 

^. <?. Reparametrization of a branch corresponding to 9 = 1 yielding a branch for 9 < 1 

We recall that 9 G [&{p,d), 1] with ^{p,d) = d^. On the cylinder C := R x S'^'^, we 
are interested in the minimization problem 1/Kckn(6', A,p) = infi,gHi(c)\{o} Qa^]- Any 
minimizer v can be taken nonnegative without restriction and solves the Euler-Lagrange 
equation 



-Av+l 



\l-0)t[v] + A\v-l-^^vP-'=O in C, (6) 



with t[v] = -j^r^^ and Qa[^] := Q^M - ^^ 
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If 6' = 1 and A = /i, any nonnegative minimizer u solves the Euler-Lagrange 
equation 

-^u + ^iu- ^^Jj wP-^ = in C. (7) 

II"'IIlp(c) 

For any 6 E ['d{p, d),l), solutions of ^ are in one-to-one correspondance with solutions 
of (It]) through the relation u = v and (1 — 9) t[v] + A = 9 jj. Summarizing, if for any /i 
we have a solution m^ of ([7| , and if we define 

f \ +r 1 A r \ II"mIIl2(c) 
T{fi) := t[u^\ and z/(/i) := - — rr^ , 



we can parametrize A and 



''^/^llLP(C) 



\V7 l|2 J_ A II l|2 V II ||2(1-^ 

IkTP 
II"mIIlp(C) 



by 

A\fi) = 9fi-{l-9)r{fi) and J^(/i) = z/(/i) 0^ (/i + r(/i))^ (8) 

For consistency, we adopt the notation u^ for solutions of ([T]) which are non-symmetric 
in the range of symmetry breaking, and m^,=k for the branch of symmetric solutions which 
are optimal in the corresponding class. 

All computations of this section are valid for the branch of non-symmetric solutions, 
but can be done in the class of radial functions as well. With evident notations, we shall 
then denote by Uf^^^ the solution of (0), and define t^,{^), ^'*(/i), Af(/i) and Jf (/u) 



2.4- Gagliardo-Nirenherg inequalities and the corresponding asymptotic regime 

In this section we investigate the asymptotic regimes corresponding to A — > oo. Let 

S,(M^) := inf kA^u?dx + ^^Au\^dx ^ 

and recall that the optimal constant Kqn = Kgn(p, c?) in the Gagliardo-Nirenberg- 
Sobolev inequality is given by 

llvy ||2i? II ||2(l-i?) 



Kgn «eHi(K'')\{o} |§d-i| V IImIP , 2p 
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An optimization of the quotient in the expression of Sp(M'^) allows to relate this constant 
with Kgn- Indeed, if we optimize M[u\ := Jj^^ |Vup dx + yu/j^d I^P dx under the scaling 

d 

A I— > u\{x) := Xp u{Xx), we find that 






achieves its minimum at 



^ , ^/u I|w||l2 



1-^ ||Vn||L2, 
so that 

ATK] = ^-'^ (1 - ^)-('-^^ II V«l|^.V)ll«llL^i? -"'"' 



thus proving that, with the choice /x = 1, KqJ^ = 'i?'' (1 — ?9)^~''Sp(]R'^). For any /i > 0, 
if M^ is the solution of (JTl) and if it is a minimizer of 1/Kckn(1, A,p), we know from [5l 
Theorem 1.2] that as /i — )■ cxd, 

(r(/i) + /x) K/^) = Q,K] ~ S,(M^) /-^ 

If n is an optimal function for Sp(M'^), we also know from the above computations that 
Aj, = 1, that is, 

1 = A: = — ^ —— and so . - 



1 - ^ r(/i) ' /i 1 - ^ ■ 

Hence, 

i^(/x) ~ (l-^)Sp(M'^)/x-'^ as /i^oo. 

Consider now the case 9 > i!) = 'd{p,d). According to ([s]), that is, by definition of 
J^ and A^, we obtain that 

as /i — > oo. Hence we have proved the following asymptotic result. 
Proposition 1 With the previous notations, for all 6 > •& = t?(p, d), 

limA^jV) = ^(l-^r^K-. 
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Moreover, the parametric curve fi H- (A^(yu), J^(/i)) is asymptotic to the curve 
9^{l-^{p,d)f-^M 

^^ {e-^{p,d)y-^M ^'^^^^ 

for large values of fi or, equivalently, for large values of A = A^(/i). 



See Figs. [lj-[3j for some plots of the curves n i— >■ (A (/x), J (/x)) for various values of 6 
and how these curves can be compared with the ones corresponding to the asymptotic 
regime as described in Proposition [l| 

The limit case 9 = d = 'd{p^d) is of particular interest. Indeed, according to [16], 
Gaghardo-Nirenberg inequalities play a special role. See Fig. |6| First of all, since 
I/KcknI^?, A^(/i),p) < J'^(/i) and using the fact that A i— )■ Kckn(^,A,p) is a non- 
increasing function of A, we recover the known result that 

Kgn<KcknWp,c/),A,p) VA>0. 

Such an inequality has deep implications on the existence of an optimal function (see [TB] 
and in particular p^ Theorem 1.4]): either the inequality is strict and there exists a non- 
trivial optimal function for ([I]), or there is equality and a non-trivial optimal function 
may exist only if A = inf{A > : Kqn > Kckn(^(p, d), A,p)}, but certainly not for any 
larger value of A, if the above infimum is finite. 

In our setting, we can define /xgn := inf{/i > : J^{^j) < Kqn}, with {} = {}{p, d). 
It is granted that /icN > 0. Either /^gn = oo and there is always a minimizer, or 

KcKN(^(p,rf),A''(/x),p) = J\^l) V/X G (0,/iGN] 

and there exists a non-trivial optimal function for ([T| if /i < /^gnj while 

KcKN(^(p,c^),A''(/i),p) = Kgn V/i G [/iGN,00) 
and there is no optimal function for ([I]) if yU > /xgn- 

3. A formal expansion at the bifurcation point 

In this section, we determine the behavior of the branch of non-symmetric positive 
solutions that bifurcates from the branch of the symmetric ones in a neighborhood of 
the first bifurcation point. Consider first the case 6 = 1 and denote by m^,* the positive 
symmetric solution of 

-An + /iM = M^"\ (9) 
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so that Q^[u^^^] = \\u^i,*\\i^(c) ~ I/^cknII'/^'P)- Notice that if m is a solution to (9), we 



still have Qfj,[u] = \\u\\^pi^s even if u is not symmetric. 



M,*llLP(C) 
||P-2 
ML^-J — II"'IIlp(C) 

We will search for minimizers of Q^ in a restricted class of functions depending only 



on the variable s (see Section 2.1) along the axis of the cylinder and on the azimuthal 
angle ( of the sphere because of the result on Schwarz foliated symmetry of [17]. This 
guarantees that we are in the right class for minimizers when 6 = 1. For ^ < 1, no such 
result has been established in the literature but we will work in the same framework. It 
is indeed straightforward to cheek that the same result holds. 

3. 1 . Expansion of Q^ at order two 

Let /i be the first non-constant spherical harmonic function, i.e. the eigenfunction of 
the Laplace-Beltrami operator on the sphere S*^"^ corresponding to the eigenvalue d—1 
and denote by /2 the next one (among the ones depending only on the azimuthal angle 



C), with corresponding eigenvalue equal to 2d. See Appendix A. 4 for details. In the 



limiting regime corresponding to £ — > 0, an expansion at order two in e gives 

D U 1 lU \\P ^^ ' II l|2p -T uyt ), 

^fi[Uij,,*\ 1 1 "M,* 1 1 LP (C) II'"a«.*IIlp(C) 

where g[/i, ip] := L (| Vy^p + /i |(y9p — {p — 1) u^~^ Iv^P) dy. By minimizing the term of 
order two, we find that 

— 7S~^ ] -*- ~ ^ — n \\p as £ -^ U , 

where Hf^ := —^+fi + d— 1 — {p — ^)u^~^^ is a Poschl- Teller operator whose 
lowest eigenvalue is given by Ai(/i) = d— 1 + fi— ^fxp"^, and such that ip = ipi fi 



is the corresponding eigenfunction (see Appendix A. 2 for details). To fix notations, we 
normalize (p by the condition 

II ||2 II ||2 II IIP 

IIV^IIl2(c) — II'/'iIIl2(k) — II^m,*IIlp{c) ' 
which slightly simplifies some computations below. This shows in particular that 



I IIP 
I"'m.*IIlp(c) 



Ai(/i). 



See (A. 2) for an expression of ipi. 

As in [1], let /ips be such that Ai(/iFs) = 0, that is 
. d-l 



p2 _4 
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see 



Appendix A. 2 for details). For any /i > yUps, we have 

Ai(/i) = -^(p'-4)(/i-/iFs). (10) 

This determines the 0{e^) term. Now we want to investigate the behavior of Q^ in a 
neighborhood of /i = yUps and therefore need an expansion at higher order. 

3. 2. Expansion of Q^ at order four 

Our purpose is to build a formal expansion of the branch of positive solutions of ([9]) that 
bifurcates from the branch /i h-)- u^^^ at M^^pg,* =: mfs and satisfies Q^[M/i] < Q^[ii^,*]. 
For /i in a neighborhood of /xps, we look for solutions of ^ of the form 



u 



M 



u^^^+eip + r]il} + o{e) + o{ri), (11) 



with e > and rj = o{e). The fact that such an expansion can be built should be 
justified in a rigorous proof. We will not do it here and keep our computations at a 
formal level. Here (p = ipi fi has been determined above. Recall that (fi is a function 
depending on s only. For convenience, let us define 

J^ MJi,* dy 

Since we are interested in functions depending only on the azimuthal angle (, we 
indifferently use u G S'^^^ or (^ G [0, tt] with a slight abuse of notation. We consider 



the sequence {fk)ken of spherical harmonics depending only on (. See Appendix A. 4 for 
details. We denote by tpk the decomposition of ip in spherical harmonics: 

ip = k^ M^,* + y^ ^fc /fc , 



A:>0 



with ^fc(r) := / ^{r,u) fk{u:)dv{u) VrG[0, cx)), 

where dv{ijj) is the uniform probability measure on the sphere. Here we have chosen ipQ 
in such a way that J^ u^~} i/jq dy = because 

, _ Ic K~*^ ih ^1^,* + i^o) dy _ , Ic K~*^ ^0 dy 
^~ f u^ dv ~ ^ f u^ dv ' 

We know that q[n, ip] = J2k>o ?[/^' ^fc M > where 

q[f^, i^k fk] = / I V(^fc /fc) r c^y + /i 'ipldy-{p-l) / m)^7,2 Vfc 0??/ . 
Jc Jc Jc 
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With fik := fi + k {k + d — 2) , for any k > 2, we get that 



13 



g[/^, V'fc fk] = \^'k\ dy + fik i'k dy-{p-i) / u^,^ i)f^ dy 

is nonnegative for /i — /j-fs > 0, small enough, and positive unless ipk = 0. Lengthy but 
straightforward computations show that, at least at formal level, 

/"^"^J 1 = ae'^ + be'^ + ce'^r] + dr]'^ + e€r]+o{e'^ + r]'^ + \a\ e'^ + t] e'^) 



(a rigorous justification of the order of the remainder terms requires some additional 
work), where 

a(/i) = r p , = Ai(/i), 
j^ u^,* dy 



p — 1 



c{ix) 



- Al(/^) r p , + T (P - 1) (P - 2) 

J^ W^^=i, fly 4 

- TTT (P - 2) (p - 3) .^p — , 
12 j^ mJ;,* dy 



IcK,*^'P^dy 

Jc K,* dy . 



L M^,=^ rfy 



- (p - 1) (p - 2) 



IcK,*^'P^'^dy 

Jc K,* dy 



d{ji) 



L M^,* dy 



+ (p - 2) fe 



i/) ) 



e(/i) = 2Ai(/i) 



!c^'^ dy 
J^ u^,^^ dy 



With no restriction, we may require that (p is optimal in the direction /i, that is 



(fip dy = . 



(12) 



Such an assumption is of course formal and a complete justification requires delicate 
estimates. This amounts to require that e(/i) = for any yU > /ips- 



According to (10), we get 
1 



fl(/i) 



(P -4)(/i-/iFs) 
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Using J^lfil^du 



d+2 



see 



Appendix A.4), we obtain 



4 6(/i 



FSj 



ip - 1) (p - 2) 



(p-1) 



JM 



U 



P-2 ,„2 
FS 



Lfl ds 



m 



tl-pQ OjO 



d{p-3) 4^3 ¥^tds 



All above integrals are computed in Appendix A. 3 and allow to express 6(/iFs) as 

4(rf-l)2p3(p_i)2[2p(5p-6) -rf(p2-i6p+12)] 



6(/x 



FSJ 



(d + 2) (p + 2)2 (p - 2) (3p - 2)2 (5p - 6) 

As for the terms which depend on rj, we observe that they sum as 
e^ 7] c{fi) + rj^ diyjj) + er] e{fi) 



q[fx,ip] 



\u 



+ ip-2)kl 



■m,*IIlp(c) 



+ e^ri{p-l){p-2) 



- 2 Xi{fi) k^ e"^ 7] 



.1 lip l^^ 



\u 



■M,*llLP(C) 



\u 



IP 
a«.*iIlp(c) 



Using the fact that /i = /o + i^{d) J2 (see Appendix A.4), it is straightforward to observe 
that the minimizer is given by 



^ = k^ M^,* /o + ^"0/0 + ^\ /l + ^2 /2 

while 'i/'fc = for any k > 2 and solves the Euler-Lagrange equation 



(13) 



+ 



2t] 



p^p-l)ip-2) 
3p-2 



/i-2Ai(/i) <7;/o 



- (p - 1) (p - 2) <-=^ ^? (/o + K(,) h) 



+ £^1/1=0. 



Here we used the fact that /^ ^^,*^ V^^ '^2//ll'", 



't^,*\\LP{C) 



3p-2 



see 



Appendix A. 3 



for 



details). Constraint (12) is taken into account through the Lagrange multiplier £, and 



a complete justification should take into account correction terms of higher order: see 
Remark |2l 



The three components ipo, i'l and 4'2 satisfy the equations 

/3Mo 



2 [cosh(ps)J'' 



/3^5o 



[cosh(/3s)]^p-2 



[cosh{f3s)f p-2 



2[cosh(/3s)]^ 

,// , P(P~1)/^ , 

2 [cosh(p s)J2 



/3^i?2 



[cosh(/3s)]^ p-2 



Symmetry breaking in PDEs 15 

with yUi=yU + (i— 1, /i2 = /i + 2(i and 



2 



B„ = _(p-l)(p-2) ^ 

^2 = ^ (p - 1) (p - 2) K(rf) -^ = fi:(d) 5o , 
Recall that 0^ = \{p- 2)^ fi. 

Multiplying the equation for ipi by (pi and integrating by parts we get 

Ai(/i) / (fi^jJids + C / \(pi\'^ ds = . 



Using Assumption (12), this proves that C = 0. This implies that ipi is an eigenfunction 



of H^, with eigenvalue Ai(/i). Since Ai(/i) is simple, we find that ■i/'i = by (12). 
We may next observe that by taking 

tpo{s) = AoXo,p-i{l3s)+BoXo,2p-3i/3s) and ?/'2(s) = ^2 X2,2p-3(/3s) , (14) 

the problem is reduced to the set of equations 

_a/' -I- ^>^°'P-i _ 2p(p-i)xo,p-i _|_ „,,p-i _ n 

A0,p-1 "T (p-2)2 (p-2)2(coshs)2 "^ "^ — U) 

_^// _|_ 4xo,2p-3 _ 2p(p-i)xo.2p-3 _ „,,2p-3 _ n n c;^ 

A0,2p-3 "T (p_2)2 (p- 2)2 (cosh s)2 "^ " ^' K^'^) 

_^ll I 4/J2X2,2p-3 _ 2p(p-l)x2,2p-3 _ ,,,2p-3 _ n 

A2,2p-3 "T ^(p-2)2 (p-2)2(coshs)2 ^ ~^' 

2 

where t(;(s) = (coshs) p-^. Since 

-^" + ^^^^ ^P(^-^)^ +^^^^-^ = 

(p-2)2 (p-2)2(coshs)2^p-2 



it follows that 

p-2 
2p 

We may notice that the equations for Xo,p-i, Xo,2p-3 and X2,2p-3 are all independent 
of k^. Moreover, since f^u^^^ ■j/'q c^-s = 0, we get 



I.e. 



/ w^ ^ {Ao Xo,p-i + Bo Xo,2p-3) ds = 0, 

JR 



n A I X0,p-1 J I D / X0,2p-3 , 

= ^0 / ^j^ ds + Bo / ^j^ ds 

(coshs) p-2 Jr (coshs) p- 



-2 
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Remark 2 The decomposition (13) ofip is formal, because it is done up to higher order 



terms. Hence the above equality only holds for /i = yU^S; ^s can be checked by computing 



Xo,p-i w"^ ^ ds 



p-2 
2p 



w^ ds 



p-2 
2p 



and, using (15), 

bo,p-i 



X0,2p-3 



wP-^ ds 



(also see 



Appendix A.l). With A 



2ri 



•Ip-i 



2p-3 J V ^ I 

Xo,p-i w P ds = - \p . 

3p — 2 

1^fi( ^^^^.7„^irV -2AiM) and Bo 



3p-2 



{p — 1) (p — 2) " „2 , we find that 



/32 



p-2 p-2 

Aq —: Ip — Bq Ir, = 



2p 



3J9-2 ^ 



holds if and only if Xi{fi) =0. As we shall see below, this is consistent with our 
expansion in terms of powers of e and rj because for ft > /ips, close enough to fips, ^lif^) 
corresponds to a term of higher order. 



The reader is invited to check that 

X0,2p-3 = 



1 p-2 , 

- \2w — w 



4 p 



As a consequence, one can compute 



P-n 



» -/ 

Jr 



Xo,2p-3W ds = - 



1 p-2 
\p-\ 



;2i2-g 



p{p-2) 



2(p-i)(p+2: 



l2.(16) 



Altogether we have found that 

e^ r]c{ii) + rf' d{[i) + er]e{fi) = r]'^ Q[ip] - e"^ r] L[ 
up to higher order terms in e, t] and (/i — /^fs)) where 



1 1 "M.* 1 1 LP (C) 



and 



ip-l)ip-2) 



Ic K,*^ ^^ dy J^ uP^/ ip^ iP dy 



kM,*r 



M:*llLP(C) 



iw/^,*r 



LP(C) 
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are respectively quadratic and linear with respect to if). Since we can multiply ip by 
any positive constant v and rj by 1/z/ simultaneously without changing the value of 
V^ QW\ ~ ^"^ V LW\-i W6 do not introduce any restriction if we make the ansatz 

V = e\ (17) 

thus making the sum of the two terms equal to £:^(Q [■?/'] — L[ip]) and ip independent 
of e, at least at leading order. Moreover, if ^/^ is a minimizer of Q[ip] — L [?/'], then it is 
straightforward to check that 2 Q[iIj] — L[iIj] = (multiply the Euler-Lagrange equation 
by ip and integrate). Altogether, we have found that 

e^ric^jj) +rj^ d{n) + erie^jj) = --e'^ L[^]<Q 

up to higher order terms in e and for /x — /xps small enough, if ^ is a minimizer of 

Q[il)\ — L ['?/'], that is, 

^^''""^^ l-\i{^)e''-h{^i)e^ = -\L[i,]e^ + o{e^). (18) 



Q,KA '^"^ ^"^ 2 



At this point, we may notice that the solution u^ has not been normalized. Multiplying it 
by a constant would not change the value of Q^^fu/^]. If we want it to be a solution of ([9]), 

■J r 1 II IIP" 2 

'-A^] — II'"IIlp(C) 



then this implies that Q^[m] = ||^i|lLP(ci ^^"-^ ^^ may therefore impose the corresponding 



constraint, i.e. 



Vm^I dy + h u^ dy = u^ dy , 
c Jc Jc 



without changing the equations written order by order (in other words, the Lagrange 
multiplier associated to this constraint is zero). Written in terms of (p and ip, at lowest 
order, that is at order e^, this constraint amounts to 

^ (l V^P + /x l^p - P^P^ u^-' l^p) dy-{p- 2) k^ ^ u^^^^ dy = 0. 

Hence by taking the limit as /x — )■ (/iFs)+ and observing that Ai(/x) = 0(/x — /xps), "we 
find that, for fi = /ips, 

k =_lr . Jc^'l^l'dy _ 2p^p-i){d-i) 

"" 2^^ ^ /c<*rfy (p-2)(p + 2)(3p-2)- 

The explicit value of k^p will however not be needed later, because of cancellations that 
occur in all subsequent computations. 

Next comes the observation that, as long as we are interested in computing L [■?/'], 
we do not even need to normalize m^ nor to compute k^. Indeed with ip := ip — k^ m^^^, = 
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\v F '"'^ 11?/ F ii- 11^ 



''M.*llLP(C) 



\u 



M.*llLP(C) 



where ■?/' is fully determined by the coefficients Aq, -Bq, and -B2, and Equation (15). This 
also determines 



(p-l)(p-2) 



/c <,*^ '/'^ V" dy 



\u,. 



>>*llLP(C) 

which is not known explicitly but is independent of k^ and can be computed as 



/c K,* ^^ ^ dy _ Ao p-2 Bo p_2 bo^ 



\u 



m,*IIlp{c) 



a 2p 



ao H a^ 

a 



p-3 



Bo 



vP-2 



a ^ ' I 



2,2p-3 



with 



ao : = 



/c K,*^ ^^ c?y p2 



/i 



/^ M^,* (i|/ 3 p - 2 



) D0,2p-3 



X0,2p-3 



tt;2p-3 c^s 



and b 



2,2p-3 



X2,2p^3w'^P ^ ds. 



Notice that here we have taken advantage of the facts that xo,p~i{(^ ■^) = ^~ „ '^{^) and 
ip^ = ifl (/o + K(d) 72) . Using the expression of Xo,2p-3, we can also compute 

1 p-2 



-'0,2p-3 



. 2 / w^^'P-^^ ds- I w^P-^ds 
{p-2)p{3p-4:) 



(j9-l)(3p-2)(5p-6) 



U. 



Altogether, with y := ^'^.^ ^ , we have found that 



4(d-l) 



2 (p-l)p^ 
(p + 2)2 



{p-2)p rf-1 p-1 



(3j9-2)2(5p-6) rf + 2 (p-2)^ 



We have not been able to find an explicit expression for b2,2p-3, but we can prove 
that this is a positive quantity and even give an upper bound. According to [18, p. 74], 
the lowest eigenvalue of the Poschl- Teller operator ~j^^Uq (coshs)^^ is given by 

1 



Ao=-Vl + 4f/o---t/o,, 

if we assume that Uq is positive. Here we have that Uq is given by 

2p(p-l) 



Uo{p) :-- 



{p - 2y 
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and the reader is invited to check that 

4/i2 



(j{p,d) := Ao + 



/xfs (p - 2)^ 



^v/rT4f^-i-f/o(p) + ^ 



2 + 



d (p2 - 4) 
d-1 



is larger than 1 for any p > 2 and any d >2. As a straightforward consequence of (15) 
we deduce that 



^(P,f^) I|X2,2p-3|Il2{R) < / X2,2p-3W^P ^ ds < 



I II II 2p— 3|| 

|X2,2p-3||l2(K) 11^ ||l2( 



and, finally, 



J2,2p-3 



< 



a{p, d) 



«;2(2p-3)rfs 



16p(p-l)(3p-4) I, 



(3p - 2) (5j9 - 6) (7p - 10) a(p, d) 



3. 3. Optimization 



Collecting the above estimates and using (18), we get 
provided /i — /ips = Oie"^). Let 



%)-2^[ 



6(/x)--LH 



-1 



Im=mfs 



so that, in a neighborhood of /x = /xfs, 



2/.K 



TtPt = ^-W-^) ( (/^-/^Fs)e^ - ;^ 



e' + oie^) 



+ o(e^) 



is optimized, up to higher order terms, by taking 

e^ = e^(/i) ~ Cp^d (yu - yUFs) as /i -)■ /iFS+ 
if Cp^rf is positive. Altogether, we have found that 



Q>, 



fi["'fn 



Q,h 



1- 



p^ 



Cp,d (/i - /iFs)^ + 0{ifl- fJ'Fsf) ■ 



(19) 



(20) 



ML^-Mr 
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All above computations are valid under the assumption that Cp^d is positive, but this 
is not a priori granted. With the estimate of b2,2p-3 that has been found at the end 



of Section 3.2, we can a posteriori give a sufficient condition for the consistency of the 



method. Since 



< L 



approx L 



4(rf-l) 



{p — 1) p^ 

Ip + W 



{p-2)p 



+ 2 



d-1 p-1 



(3p-2)2(5p-6) d + 2 {p-2) 



with y : = 

-'^approx [y , 



16p(p-l)(3p-4) 



(3p-2) (5p-6) (7p-10) (T(p,d) ' 

< 2b{fips)- Moreover, we have 



cp,.<c-r-=-(p^-4) 



we know that Cp^d is well defined and positive if 



6(/i) 



^approx [ 



Ia»=mfs 



at least as long as i^approx[V'] < 2 6(yUFs)- This last condition holds Up is contained in a 
non empty interval (2,papprox) C (2,2*), where Papprox is determined as the largest root 
of a fourth order polynomial, namely p i— )■ | (54 — 227 c? + 103(P)p'^ — 4 (18 — 37 d + 
67rf + 46t/2)p2^ l6{d + 3){5d-3)p-2A0d{d+ 1). In practice 



25 d 



'2^3 + 1(63 



d !-)• Papprox('^) is closc to 2d/{d — 2) and converges to 2 as (i — ;■ cxo. See Fig. 11 



3.5. Expansion of T{fi) around ^-ps 

Let r(;.) = t^] = ^-^^^^ and 

p-2 



(21) 



We can notice that T(/iFs) = /3^ J2/I2 (see Appendix A.l and Appendix A. 2) so that 

p-2 4(rf-l) 



T{fi 



FSJ 



yUps 



p + 2"^^^ (p + 2)2 ■ 
With the above expressions in hand, we can now compute the derivative 



- — (fi — /ips)- By expanding the expression 



From (21) we know that ^['U^,,,,] — t[uFs] — ,2 

T{fi) - r(/iFs) = t[Ufi] - t[ups] = t[u^^^] - t[uFs] + t[u^l] - t[Uf,^*] in powers of e with 



u^ = u^^^: + Eipi fi+ e^Tp + 0(5^) where </?, ip have been chosen in Section 3.2 , we get, up 
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(Ai(/i)-/.-tKJ)y^j^^— + (p - 1) 

jc iJ-,* y Jc fj-,* y 



IcK,*^y Ic^l,*^y 



and, by computing as above, we find that 



r'(/i 



FSJ 



p-2 
p + 2 



+ c 



p,d 



c ^1 -^ + (p _ 1) ic <s' ^? c?2/ 



p + 2 J^ulsdy ' ^"^ ^^ Ic'^Fsdy 
4pyUFS /cWfsV'o c^I/ 



p + 2 /^ m2 ^ rfy 

because we notice that the terms involving k^ canceL Hence, using (13) and (14), we 
have found that 






4pVfs ^ 2(p-1)pVIs 



(p + 2)2 (p + 2)(3p-2) 

2(p-2) Ao 4p/iFs5obo,i 
p + 2 a p + 2al2 



p-2 
p + 2 



vA 



8p{d-l) 
{p-2){p + 2y 



1 2 



(p-l)p(p + 2) bo,i 
2(p-2) I2 



Here the coefficient bo,i is given by bo,i := J^Xo,2p-3'UJ ds. Using (16), we obtain that 
r'(/iFs) = ^—: + , \,w , ' Cp,rf . (22) 



3.6. Expansion of ui^fj,) around fj,Fs 



\u,. 



Let us consider z/(/i) := ||^^||l2(c)/ ii"'miilp(c) 

I/(/i) - I/(/XFs) = 



Again we can write 



2 II 112 

^m,*IIl2(c) II^fs|Il2(c) \ , / / n II^m,*iil2(c) 

' + z/(/i) 



M 



M:*llLP(C) 



FfsI 



LP(C) ^ 



M 



With /3 = ^2^ y^, using expressions that can be found in Appendix A.l 



M,*llLP(C) 

we see that 



, . II'"m,*IIl2(c) 0-2=^ I2 



\u 



■m.*IIlp(c) 



2/p 



Kp fi 



P-2 
2p 



(23) 



with.,= (^)^(|^^ng(. 



i^*(/^ 



FSj 



and hence 



I l|2 II ||2 

Fm,*IIl2(c) IFfs|Il2(c) 



hm . 

/i-s-MFS /i — /ips \ ll^/i,*llLP(C) IpFsIIlp(C) 



p-2 Z/=,(yUFs) 



2p /ips 
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If M^ = u^^* + e({) + e Tp + o{e ), where y?, ijj have been chosen in Section 3.2, after a 
Taylor expansion we find that 



i + .mJ^^-Ip-D^'^ "'""'* 






^^^2!c^i^,*^dy ^^2 



/^ u^,* (iy 



Again we may notice that the terms involving k^ cancel and, based on (13) and (14) 



we arrive at 



v\^^ 



FSJ 



P 



Z^*(/i 



FSJ 



2p/iFS 



+ C„d 



P/Ups 



p{p-l) 



+ 2 



p-2 
2pnFs 



+ c 



'P,a 



p + 2 3p-2 

Bo bo,i Bq bo,p-i 
a I2 a \p 
2 p{p-l) 



P/^FS 

+ 2 



p+2 3p-2 
Bo /bo,i p-2 
a V I2 3p-2 



Claim 3 ili ^/ie bifurcation point fi = /ips '"^e g'ei i/ie following inequalities. 



Ki^^ 



FSj 



+ 



r:(/i 



FSj 



^'(/i 



FSj 



+ 



^'(/Ufs) 



l^*{l^Fs) /^FS + T"*(/^FS) '^*{l^Fs) fJ'FS + nifJ'FS 



0. 



Proo/. Recall that Q;^[m^,*] = z/,,(/i) (r*(yu) + /i) and Q^[m^] = z^(/i) (T(/i) + /i), so that 



d 



QfiFsi'^Fs] dfi 



'^fj,["'ii,*\ 



Kil^Fs) , l+^i(AtFs) 



according to (21) and (|23|), and 
1 



d 



Q/^fsNfs] C?/i 



QmK 



^''Fs z/^ (/ips ) /ips + r* (/ips ) 



'^'(/iFs) 1 + r'ifiFs) 



=ws z/(/ips) /ips + r^fiFs) ' 



According to (20), these two quantities are equal, thus proving the result. Alternatively, 



the identity can be proved directly using the expressions of r' and z/' established in 
Sections 13.51 and 13.61 D 



3. 7. Reparametrization of the branch for 6 < 1 

Now we are in position to study the local behavior of the branch of the solutions to ^ 
parametrized by /i close to the bifurcation point, that is, for /i in a neighborhood of /ips. 
More precisely, we are interested in the monotonicity of /i h> A^(yu) and the behavior 



Symmetry breaking in PDEs 



23 



of yU t-T- (A (/i), J (fi)) in a neighborhood of /i = /ips- According to ([8j), we know that 
A%fi) = 6fi- {1- 6) T{fi), so that 



{A'y = e-{i 



T 



can be computed at /i = yUps using the expression of r'(yUFs) that has been computed in 
Section I3.5[ Hence we find that 

(A^y(/iFs) = ^(l + r'(/iFs))-r(/iFs). 

Claim 4 //r'(/iFs) ^s positive, then we have that ^ A^(/iFs) < if and only 

^'(/iFs) 



e<Mp,d):-- 



1 + ^'(/iFs) 



Notice that with this definition, 'i92(p, d) is defined for any p G (2, 2*) and any d> 2. As 
long as 'd(p, d) < '&2{p^ d), (A^)'(/iFs) is negative if 6* G {^{p, d),i!}2{p, d)). In all numerical 
examples that are under consideration in this paper, we find that r'(/iFs) is positive. 



This is of course automatically the case if Cp^d itself is positive, because of (22). 

We recall that J^fi) := 6'" (/i + r(yu))^ u{ij,) and A^(/i) = 9 fi - {1 ~ 9) t{^i). Hence 



{log j'y 



V 



fi + T 

and a similar formula holds for Jf. At // = /xps, we can use Claim [s] and get 

^(l + r'(/XFs))-r'(/iFs) (A7(/^Fs) 



(logr)'(/iFs) 



/ips + T"(AtFs) 



/iFS + T"(/iFs) ' 



Let us define 



iJi' (j:) 



e\/ 



■ (A7 (AD'- 

Claim 5 Assuming that remainder terms are all of lower order in the above compu- 
tations, we have 



S\ix 



FSJ 



0. 



In other words, the parametric curves /x i— )■ (A^(/x), J^{fJ.)) and fi h-> (A^(yu), Jf (/x)) are 
tangent at /z = /xfs- It remains to decide the relative position of the two branches in a 
neighborhood of the bifurcation point. In order to do so, let us define the function 



e-.-- 



((A 



6l\/\2 



J' 



d 



/(J7 



1 d 



f M' 



L(A7rf/iV(A7; (AD'rf/iV(AD' 
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which, up to a positive factor, is the difference of the second parametric derivatives of 
J^ and Jl. 

Claim 6 Assume that i)2{p,d) > 'd{p,d). The parametric curve fi i— ?■ (A^(/i), J^(/i)) is 
concave and below the parametric curve fi t-)- (A^(/i), J;f(/i)) for ji — /ips > 0, small, 
if 9 > 'd2{p,d) and ^^(/Ups) < 0. The parametric curve ji t-)- (A^(yu), J^(/i)) is above 
the parametric curve /i i— ;■ (A^(yu), J^(yu)) for jj, — /xps > 0, small, if 9 < 'd2{p,d) and 

e'(/iFs) < 0. 

Now let us discuss the range of positivity of ^^. An elementary computation shows that 

{jy {j'y {Ay {jy f{A'y' ^ 



9) T{fi) and {}2{p,d), we 
'-^^^(P^. From (S and 



^ J^ (A^)' J^ J^ \{Al)\ 

because (A^)" = 0. By definition of A%ij) = 9 jj - {1 
get that (A^)"(^ps) = -(1 - ^)r"(/xps) and (A^)'(/ips) = 
Af(/i) = 9fi- {1 - ^)r,(/i), we get that (Af)' 



2pe-(p-2) 



p-2 



p+2 



From (21) and (23), 



{JtY 



2pe-{p-2) 



we know also that Jf(/i) = /tp {2p9/{p + 2)) fi '^p , so that ^ — 2 
(log jy' = — ^2~ 2 • As a consequence of Claim pi we can write that 



and 



id'Yif^ 



j^(/. 



Ps) ^ (J^Vfs) (A^IVfs) 
Ps) J^XfiFs) (A^)'(/^Fs) 



According to (20), we have the identity 

J' 



1 



{p' - 4) Cp,d = log 



Jf 



(/xps) + (1 - ^) lo, 






|A'=MFS 



which allows us to compute 

(^^)"(/iFs) 
J^(/XPS) 



- - (p2 - 4) cp,d + 






-{1-9) 


P + 2 „, ^ / P + 2 ,^ ,, \^ 1 






-(1-^) 


2p/ips V2pyUFS / yUps 




of Claim |5| we also have 


(^')'(/^Fs) ( 


A')"(/^FS) (^i')'(/iFs) (A^)"(/^FS) 





(A^)'(/ips) J^(/ips) 



Jffe) (Af)'(/iFs) 
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p + 2 



{l-e)r"{^^^s) 



Collecting the above identities, we can compute the value of ^^(^ips) as 



eif^ 



FSJ 






4) Cp^d 

+ (iogj:)"(/iFs) 

1-e 



(A7(/fFs)Y 



/^FS 



1- 



p 



2^^(i + ^W)))' 



The cancellation of the terms involving r"(yUFs) is a remarkable fact. By definition of 
■&2{p,d), we get 

e'(/iFs)= -7(p'-4)cp,, 



2pd-{p-2) 



2p//|s 




1 



p + 2 



A'-FS 



p + 2 1 

2p l-i?2(p,o?) 



2pe-{p-2) l-^2ip,d) 

2 



and finally arrive at 

e'(/iFs) 






4)c P + 2 (l-g)(2p^2(p,rf)-(p-2))^ 

^''''' 4pV^s {l~MP,d)f{2pO-ip-2)) 



At this point, we can observe that 'd(p,d) > ^— . The reader is then invited to check 
that the function 6 h-j- ^^(/ips) is nonin creasing on [^{p, d), 1] and 



e 



MP,d) 



(/i 



FSj 



-{p - 4) Cp,d + — ^^^2 



4pVfs 



^2(P,C^) 



because of (22). Recall that the positivity of Cp^d is required in (19). 

Claim 7 If'd2{p,d) < 'd{p,d), ^^(/ips) is negative for any 9 G {{}{p,d), 1]. Otherwise, if 
Cp,d is positive, ^^(yUps) is positive if 9 G ['!?(p, rf), '(92(p, c?)) and negative if 9 G {^2{p,d), 1]. 



Summarizing, the expansion (11) and the subsequent computations are valid, and 
make sense, as soon as the coefficient Cp^d, whose expression is established in Section 3.3 
is positive. Then for any 9 G {m.ax{'d{p,d),'d2{p,d)},l), the curve of the energies of 
the non-symmetric solutions of the Euler-Lagrange equations is concave, nondecreasing 
as a function of A in a neighborhood of the bifurcation point, and below the energies 
of the symmetric functions. If '&2{p, d) > d{p, d), then the curve of the energies of the 
non-symmetric solutions is above the energies of symmetric functions in a neighborhood 
of the bifurcation point if ^ G ['d{p,d),'d2{p,d)). To decide whether Cp^d is positive or 
not relies on numerical computations. However, the estimate of Section 3^ shows that 
this occurs at least in a large subinterval of (2, 2*). 
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4-1- Symmetric and non-symmetric branches, and their asymptotic behavior 



In [12] the branches of solutions which bifurcate from the branches of symmetric 
solutions at their smallest possible bifurcation point have been computed numerically. 
For completeness, we start by exposing some of these numerical results, which are the 
main motivation of the present paper. The branch of symmetric solutions is explicit. 
The branch ft h-)- (A^(/i), J^(/i)) bifurcates from the symmetric ones at /x = yUps and 
is computed numerically. The algorithm is based on descent techniques and on an 
iteration scheme which allows us to compute the branches of solutions by continuation. 
We carried out the computations for dimension d = 5 and various values of p and 6. 
We have of course no guarantee that the solutions that we have computed are the 
optimal ones, but at least the values that we have found are fully compatible with what 
is expected. In particular, the curve of the computed estimates of the best constant 
is a concave increasing function of A, which can reasonably be expected to coincide 
with A H- Kckn(^) A,]9). Moreover, when 6 approaches d{p,d) from above, the curve 
A i-T- 1/Kckn(^, A,p) approaches A \-^ max{K^j^j^(6', A,p), Kqn}- Last but not least, the 
asymptotics predicted in Proposition [I] are not only correct (dotted lines in Figs flfp]) 
but provide a good upper estimate of the curve in the whole range A > 0. 




asymptotic 
symmetric 
non-symmetric 



■^bifurcation 



Figure 1. Parametric plot of ^ i-^ (A^(/^), J*(/x)) for p = 2.8, d = b, 6 = 1. Non- 
symmetric solutions bifurcate from symmetric ones at a bifurcation point /i == /ips 
computed by V. Felli and M. Schneider. The branch behaves for large values of A as 
predicted by F. Catrina and Z.-Q. Wang. 



4-2. Two scenarii 



The branch of symmetric minimal solutions, (A, m^), is known explicitly and is monotone 
in energy, that is, the function A h-)> Q^(m^) is monotone increasing in (0, +oo). In the 
computations described in [I^ we observe that the branch of non-symmetric solutions. 
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IShrS 



:j«(m) 




^bifurcation 



asymptotic 

symmetric 

non-symmetric 



a«(m) 



Figure 2. Parametric plot of ^ i-> {A^ (^i) , J^ (fi)) for p = 2.8, d = 5, 9 ^ 0.8. 
The behavior is similar to the case 6 — 1 wp to the reparametrization described in 



Section 2.3 while the asymptotic behavior of the branch for large values of A is in 
agreement with the results of Section |2.4[ 




asymptotic 
symmetric 
non-symmetric 



a'^(m) 



Figure 3. Parametric plot oi ^i ^ {A%fi) , j'^ {^i)) for p = 2.8, d ^ 5, 6 = 0.72 > 
t?(2.8, 5) ~ 0.7143. The asymptotic behavior of the branch for large values of A is in 
agreement with the results of Section [2]4j We may however notice that the bifurcation 
seems not to occur where expected. More detailed computation in the Section [43| will 
shed light on this phenomenon. 



(A^(/i), Qfi[u^]), is monotone for some values of 6 (for instance for 6 = 1), but not 
always. More concretely, for certain values of p and d, the numerical results show 
that there exists an exponent t?i = 'di{j>,d) e (^^(p, d), 1) such that for any 6 & [^i, 1] 
the branch (A^(/i), Q^[m^]) is monotone increasing. But when 'd{p,d) < 9 < 'di{p,d), 
a dramatic change occurs: see Figs. [5] and |6| For the values of p and d that have 
been considered numerically in those figures, the branch is not monotone anymore for 
/i > /ips, thus producing non-symmetric solutions and candidates for optimal functions 
in the inequalities for values of A < A^(/iFs)- This phenomenon provides an explanation 
for the results proved in [B] using rigorous a priori estimates. 



The limiting case 9 = 'd(p, d) is very interesting: see Figs. M and ^ Let us define 
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■^gnIPj d) ^^d /iGN such that 

KGN = K*cKN(%,c?),A^N(P,rf),P) and Af^''^)(/iGN) = Agn(p, ^) • 

Whether Kgn is larger or smaller than K*Q^^(J}{j),d),K-ps{j)i'^{'Pid)),p) determines, at 
least in the framework of our computations, whether 9 is smaller than ■di{p,d) or not. 
This has been observed in [13] and theoretical consequences have been established in [3] , 
in the limit regime p — ?■ 2. Before going further, let us observe that i)i{p,d) is an 
exponent associated to a global property of the branch. 

Based on our numerical computations, we are now in position to formulate the 
following alternative. 

Scenario 1. If Kgn < Kqj^jv^('(9(p, d), Afs(p, ''^(p, (i)),p), then for all 6 > ^{p,d), the 
optimal functions are symmetric for any A G {0,Aps{p,0)] and the branch of non- 
symmetric solutions is optimal for any A > Aps{p, 0). Such solutions exist for arbitrarily 
large values of A if 6* > ^(p, d): see Figs. [7]and|8| but may exist only for a finite range 
of A if e = ^(p,d) and Kgn < K^kn(^(p,c?), AFs(p,^(p,rf)),p): see Fig. g 



Scenario 2. If Kgn > ^cKN{.'&{Pid),K-ps{p,'&{p^d)),p), there exists 'di = 'di{p,d) G 
(i9(p, d), 1) such that for any 6 G [di, 1] the branch is monotone increasing. We further 



observe numerically that -^i = '§2 (see Fig. 12), where ^2 has been defined in Section 3.7 



for any 6 G [^{p, d) ,^2{p-i d)) ^ we know that the branch of non-symmetric functions is 
decreasing in a neighborhood of the bifurcation point, but also has a larger energy than 
the symmetric solutions of the Euler-Lagrange equation (for the same value of A) . Hence, 
in the critical case 6 = ■(9(p, d) we have Aq^{p, d) < Afs(p, "^{p, d)) and for any A in the 
range (0,Aqjv^] := {A,,,^^' '{^) : /i G (0,/iGN]}, the optimal functions are symmetric 
and KcKN{'^{p,d),A,p) = KQ^^('d{p,d),A,p). See Fig. [6J Moreover, in the language 
of the concentration-compactness method, according to ^E\, for any A > AQ^{p,d) the 
optimal constant is determined by the problem at infinity and Kckn(''9(p, d),A, p) = Kgn- 
From the viewpoint of the reparametrized branch, we numerically observe that both 
/i I— 7- A'^''P'^\ix) and /i t-)- J^^P''^\ix) are decreasing for /x > fips, at least for the values 
of p for which computations have been done. In the subcritical case corresponding to 
t?(]9, d) < 6 < '&2{p, d), the reparametrized branch (A^, J^) is not monotone. Numerically 
we observe that it is monotone for 9 > ^2{,P, d), hence suggesting that ^i{p, d) = ^2{p, d) 



(see Fig. 12) 



Altogether, based on our numerical observations, what decides between Scenario 1 
and Scenario 2 is the relative value of Kgn and K^j^j^('i9(p, d), Afs(p, ■(9(p, d)),p). Equality 
of these two optimal constants determines a value of p = p*{d). Numerically we find 
that j9*(5) ^ 3.001 and ■(9(p*(5),5) = ^i(p'^(5),5) ^ 0.834. For p G [p'{d),2*), only 
Scenario 1 occurs (numerical observation). For p G (2,p*(t/)) we have ^{p^d) > di{p,d) 
and Scenario 2 occurs. More precisely, the fact that the branch cannot be globally 



monotone increasing if 6' < {}2{p,d) is a consequence of Section 3.7 while the fact that 
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the branch is monotone increasing if ^ > ■&2{p,d) is a numerical observation. These 
results are fully consistent with the ones of [3j and [13]. Now let us give some details. 

4-3. Bifurcations and qualitative dependence in 9 



In Fig. |3| a careful inspection shows that the symmetric and the non-symmetric branches 
of solutions differ for values of A strictly less than A^(/iFs)- This is not the case for 9 
close enough to 1: see Fig. Ill but very clear on Fig. [SJ When 9 approaches •d{p, d), the 




Figure 4. Parametric plot of ^ H^ (A^(/x), J''(^)) for p = 2.8, d = 5, = 0.95 
close to the bifurcation point. For A < A^(/iFs), a-H solutions are symmetric, while for 
A > A^(/iFs)j non-symmetric solutions provide a better constant in the interpolation 
inequalities, exactly as for the case 9 = 1. 

branch (locally) converges to its limit: see Fig. ^ Figs. ^ and ^correspond to Scenario 2 
(for 'd <9 < -i^i), that is, to the case 



Kgn > K*ckn(^(p,c?),Afs(p,^(p,c?)),p). 

In other words, this means that p < p*{d). 

The case p > p*{d), i.e., Scenario 1, also occurs and corresponding plots are shown 
in Figs. IrHgj There we take d = 5, p = 3.15 > p*{5) ^ 3.001. 
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Figure 5. Parametric plot of /i i~> (A^(^), J^{^)) for p = 2.8, d = 5, 6* = 0.72 close to 
the bifurcation point. Non-symmetric solutions exist for A < A^(/zfs)- There exists a 
value Aq^ < A^(/iFs) such that optimal functions are symmetric for any A S (0, Aqj^) 
and are non-symmetric for A > Aq^. When A — Aqj^, symmetric and non-symmetric 
optimal functions co-exist. 
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Figure 6. Critical case 9 ~ d{p,d). Parametric plot of /x H> (A^(/j,), J^(/i)) 
for p = 2.8, d = b, e ^ ^(2.8, 5) « 0.7143 close to the bifurcation point. Non- 
symmetric solutions exist for A < A^(/iFs) but, at least for the chosen values of p 
and d, are never optimal functions for the interpolation inequalities. There exists a 
value Aqj,^ < A^(/iFs) such that optimal functions exists and are symmetric for any 
A e (0,Aqi^] and do not exist for A > A*q^. Moreover, A}.^ = A^(/iGN) with the 



notations of Section 2.4 and KcknI^^Ip, d), A,p) = Kqn for any A > Aqj,^. 
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./"(/') 




Figure 7. Parametric plot of /i ^ {^^{fJ-), J'^if^)) ioi p = 3.15, d = 5, 6 = I. The 
behavior of the non-symmetric branch is similar to the one found in Fig. IT] 




Figure 8. Parametric plot oi ^ ^ (A^(m), J\f^)) for p = 3.15, d = 5, 9 = 0.95. The 
behavior of the non-symmetric branch is still similar to the one found in Fig. [4] 



./"(m) 




A»(rt 



Figure 9. Parametric plot of /i n- (A^(/i), J^(/i)) for p 

e = 19(3.15, 5) w 0.9127, that is for the critical value of 9. 



3.15, d 
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In Section [3] we proved that the symmetric and the non-symmetric branches of 
solutions are always tangent at Aps- What happens in a neighborhood of the bifurcation 
point is therefore difficult to decide in view of the plots of the branches, especially when 
^ is in a neighborhood of ^i{p,d). To illustrate this difficulty, we may for instance 
observe that figure Fig. [5] is an enlargement of Fig. [3j Hence, in Scenario 1, we have to 



discard the possibility of other scenarii than the ones described in Section |4.2| at least 
in a neighborhood of the bifurcation point. 

The computations of Section |3] are formal for two reasons. We assume that the 



solutions in the non-symmetric branch correspond to the ansatz given by (11) and (17) 



up to various higher order terms, and we also require that Cp^d is positive. This last 



condition has been discussed in Section 3.4, but can be checked numerically within the 



ansatz: we know that Cp^d is positive and finite as long as Cp^d'^"^ i^ positive (see Fig. 10 



and Fig. fTTl for a discussion of the sign of l, ^ 
always positive. 



approxN 



, and numerically we find that Cp^d is 




Figure 10. Computation of Cp^d with d = 5 as a function of p. We observe that the 
numerical solution is positive for any p e (2, 10/3) where 10/3 is the critical exponent 



(corresponding to the plain vertical line). The estimate of Section 3.4 corresponds to 



the dotted line and holds for p < 3.2323 (corresponding to the dotted vertical line). 



The plain curves are the zeros of Cp^'J"^™ as d varies in the interval (2, 10), 



Figure 11 

except p — 2 which is a singularity. The dotted curve is given by d 



2d/{d-2) =2*. 



For a given d, 



approx 



is therefore positive for a large subinterval in p of (2, 2*). 
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Under the above precautions, we know from Section [3] that there exists a number 
d2{Pid) such that the behavior of the branch in a neighborhood of the bifurcation 
point /i = /xfs discriminates between two regimes corresponding to 6^ > 'd2{p,d) and 
9 < '(92 (p, d). When 9 < t?2(p, d), we have (A^)'(/iFs < and the contrary happens when 
9 > d2{p-,d). So, locally, the reparametrized branch is monotone increasing in a right 
neighborhood of /xps if and only if 6' > ?92(p, d). Since global monotonicity implies local 
monotonicity near the bifurcation point, if the numerical computations of the branches 
are consistent with the formal study of the bifurcation carried out in Section |3} then we 
should have that ^i{p,d) = ■t}2{p,d). It is not easy to establish a qualitative property 
such as the monotonicity, but at least we observe in Fig. 12 that for 9 = d{p,d) the 



range in p for which ■&2ip,d) > i){p,d) corresponds to the range in p for which the 
Gagliardo-Nirenberg constant compares well with the energy at the bifurcation point. 



0.15 k 



0.10 ^ 



0.05 



-0.05 - 



-0.10 




Figure 12. Comparison of tiie local and asymptotic criteria in the critical 

case 9 = i9(p, d) when d — 5. The dotted curve corresponds to the function 
p I—)- 1/Kgn(p, f^) — J^''^''^Hl^Fs), that is, the difference of the asymptotic energy of 
the branch and the energy at the bifurcation point: when it is negative, this means 
that 'di{p,d) is defined and larger than d{p,d), so that Scenario 2 takes place. When 
it is positive, this means that Scenario 1 can be expected. The exponent ^2(p, d) can 
be defined for any p G (2, 2*). The plain curve represents p i->- 5 {^2{p, d) — d{p, d)) and 
positivity indicates that, at least locally around the bifurcation point. Scenario 2 takes 
place. Hence the local (around the bifurcation point) and asymptotic (as A -^ +oo) 
criteria coincide. 



5. Concluding remarks 



In this paper we have established the asymptotic behavior of the branches for all 
9 G {i){p,d),l]. There is a good agreement between this behavior and that of the 
numerical branches, which reinforces the conjecture that the computed branches contain 
the extremals for the Caffarelli-Kohn-Nirenberg inequalities. 

We have also studied the precise behavior of the branches of non-symmetric 
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solutions near the first possible bifurcation point on the symmetric branch within the 
framework of a particular ansatz. By doing so, we have obtained the existence of a 
critical exponent "(92 (p, d) above which the branch is monotone, increasing and potentially 
optimal, and below which it is certainly not optimal in a neighborhood of the bifurcation 
point. 

The final result of this paper is the comparison of the above criterion based on 
the local behavior of the branch near the bifurcation point and the criterion based on 
the asymptotic energy of the branch in the critical case 9 = i}{p,d), using numerical 
methods. They coincide, which gives solid grounds to the alternative that has been 
numerically observed: 

Scenario 1. The non-symmetric branch is monotone increasing for any 9 G \d{p^ d),V\. 
Scenario 2. The non-symmetric branch is monotone increasing for any 9 G {'i)2{p,d), 1] 
but it is not optimal near the bifurcation point if ^ G I'&ip, d),i)2{p, d)). 

This suggests that no other scenario can take place, consistently with our numerical 
computations. Hence we arrive at the conclusion that Scenario 1 takes place when 
^2{p,d) < d{p,d) and Scenario 2 holds if d2{Pid) > d{p,d). The branches of solutions 
we have computed are likely to be optimal for the Caffarelli-Kohn-Nirenberg inequalities. 



Appendix A. Some useful quantities 

Appendix A.l. Computing integrals 

We recall that /(g) := j^ (^cotL)i ^^^ t)e explicitly computed: /(g) = -j^i+i^- An 
integration by parts shows that /(g + 2) = -tiIXq)- The following formulae are 

reproduced with no change from [B] (also see [HJ). As in [6j, with w{s) = (coshs)^p^, 
we can define 



\q ■= / \w{s)\'^ ds and ^2 '■= / \w'{s)\'^ds. 

Jr Jr 



Using the function /, we can compute I2 = /(^), \p = /(-^) = /(^ + 2) and get 
the relations 

^r(^ ) . 41, . 41, 

p+2 - 
2(p-2)' 



'' r(^) ' '^ p + 2' ^' (p + 2)(p-2)- 



As a special case, we have 



L := / (coshs)~p^ ds = ) ^—^ 
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Appendix A. 2. Symmetric extremals and linearization 

2 

Consider w(s) = (coshs) p-^, which is the unique positive solution of 

-{p - 2f w" + Aw - 2pw'P-^ =0 
on M, up to translations. The function u{s) := aw{[3 s) solves 

With (3 = 2=^ ^JJl and a = (| /x)p^, m = m^^,,, is given by 

'p \^ 



«M,*(^) = (|/^) 



1 - ■ • - 2 



cosh ( ^/i s 



P-2 



and solves Q. 

Next we are interested in computing the ground state energy of the Poschl- Teller 
operator "H^ = —-^ + d—l + ii — {p~l) m]^^^, that is the lowest eigenvalue Ai(/x) in 
the eigenvalue problem 

'H^,^l = Ai(/i)v3i . 
See [T8l [T9] for further references. The function 

Ms) := «' (cosh(/3 s))-^ = n^/2 (A.2) 

solves 

-Vi + 1 /^pVi - (p - 1) <7*^ ^1 = 

and provides a solution with 

Ai(/i) = (i — 1 + yU — - /ip^ . 
The Sturm-Liouville theory guarantees that ipi generates the ground state and 

ml T^ |T^ — Ai(_/iJ — r 112 . 

II ||2 II ||P 

II'^IIl2(c) ~ II'"m,*IIlp(c) 
Notice that the condition Ai(yUFs) = determines 

d-1 
AiFS = 4 



p2 _4 
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Ai(/i) 



4 



(/i - /iFs) • 



Other eigenvalues of l-i^ can also be computed using classical transformations and special 
functions: see [iHl p. 74]. Notice that in Section [sj we normalize the function ip = ipi fi 
in the expansion (11) by the condition ||(/?||l2(c) = ||'"a'.*IIl''(c) consistently with (A. 2). 



Appendix A. 3. Useful quantities 



Collecting results of Sections Appendix A.l and Appendix A. 2, with a^ ^ = | yU, we 



find that 



/^ uf,^^ dy 

IcK:'^'dy _f{q + 2) 



Ic ufi,* dy 



/(^) L=^ 



a 



p-2 



p-2 



3j9-2 



and 



JcK-'^Uy _ f{q + 4) 
JcK,*dy f{q) 

Ic IVy^P dy 
Ic K,* dy " I2 
1 pfi 



a 



2 {p-2) 



- 2p 
-p-2 



a 



,p-2 |p 



2p + 2 



(rf-l) + 
4(d-l) + 



P 



2p3(p_l)^2 

(3j9-2)(5p-6) ' 



p-2 

p^ (p - 2) /i 
3p-2 



/3M1- 



/(?) 



-p-2 



Tf ,, _ ,. ^ 4-V,„„ ^„ fl„J ^V,„^ Ic l^'/^P ^?/ _ 2p(p-2)(p2+p-l) 2 

it /i - /ips, tnen we nna mat ^^^2^ ^^ - (p+2)(3p-2) /^fs- 



Appendix A. 4- First spherical harmonics 

Denote by C ^ [Oj ^] the azimuthal angle and consider the Laplace-Beltrami operator C 
on the sphere S"'"^. When £ is restricted to functions on S"'"^ depending only on (, it 
takes the form 



>C/ = sin^-'^C:^(sin'^-^C'^^ 



d 

dC 



dC 



and £ is unitarily equivalent to L defined by 

lg = {l-x^)g"-{d-l)xg' xg[-1,1] 

whose eigenvalues are the Gegenbauer polynomials or ultra- spherical polynomials. The 
correspondence between the operators is simply given by 



/(C)=^(cosC) 
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and one can check that 

j^ |/(C)P sin'^-^ ^ ^^ = / 1^(^)1' (1 - ^')'^ ^^ • 
It is also not difficult to check that the first Gegenbauer polynomials are 

gQ{x) = l, gi{x) = X , g2{x) =dx^ -1 , gsi^x) = {d + 2) x^ - 3x , 
with eigenvalues respectively equal to Aq = 0, Ai = rf — 1, A2 = 2 c? and A3 = 3 (c? + 1): 

-lgo = 0, -lgi = {d-l)gi, -Ig2 = 2dg2,-lg3 = 3id+l)g3. 
On [0,7r], we consider the probability measure 

duiC) = ^ sin''-' C dC 



where 



Jo ' ' "^ r ^ 



Then 



/o(c) = 1 , /i(c) = Vrfcosc, /2(c) = y ^^^ (^ ^o^'c - 1) 

are normalized eigenfunctions in L^([0, tt], (iz/((^)): 

/'TT 

/ \fi\'diy = l V2 = 0, 1, 2, 
Jo 

with eigenvalues Ao = 0, Ai=(i — 1 and A2 = 2 d: 

-£/o = 0, -£/i = (rf-l)/i, -Cf2 = 2df2. 
We also have the useful formulae 

3d r,.,2.^ /2(rf-i) 



/^ 3(i /■'^ 

l/i|^c?z/= — -^, n(d):= j \fi? f2dv 



d + 2 
and 



fl = fo+ \l ^ ^ J2 = J0 + t^{d) 12 
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